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Abstract 

In the simplest examples of models with a discrete quark-lepton symmetry, 
an electroweak symmetry breaking sector with more than one Higgs doublet is 
necessary to obtain the correct mass relations between quarks and leptons. A 
two Higgs doublet model version has flavour- nonconserving Yukawa couplings, 
which are proportional to the masses of the quarkdepton symmetric partners 
of the fermions. We describe how flavour changing leptonic decays can occur, 
with branching ratios not far beyond that currently measurable, enabling 
investigation of the phenomenology of such models. 
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I. INTRODUCTION 



The Quark-Lepton Symmetric extension of the Standard Model involves the addition 
of leptonic colour to the symmetry group, in order that the lagrangian can be made to 
exhibit a discrete symmetry between quarks and leptons ffl]. The symmetry group then 
becomes G g e = SU(3)£xSU(3) (? xSU(2) i xU(l)x- The leptonic colour group SU(3)^ can 
be spontaneously broken at a scale as low as about a few TeV. In the simplest models, 
SU(3)£ breaks down to SU(2)', which remains unbroken, and acts to bind the exotic colour 
partners of leptons into hadron-like states. Then for the breaking of the SU(2)x, electroweak 
symmetry, the minimal choice of a single Higgs doublet as for the Standard Model gives 
rise to unsuitable relations between the masses of the quarks and leptons f2|. It is thus 
necessary to use a slightly expanded Higgs sector, and this paper will examine some of the 
phenomenological consequences of choosing the next simplest possibility, that of two Higgs 
doublets. The present paper extends the analysis given in Ref. 0. It is important to note 
however, that these unsatisfactory mass relations can be avoided by having a different form 
of quark-lepton symmetric model, such as one in which leptonic colour is broken completely, 
or by incorporating a left-right symmetry (see papers 9 and 13 in Ref. Q). 

Two Higgs doublet models (2HDMs) are usually constructed so that tree-level flavour- 
changing neutral processes are absent 0, although a generic 2HDM would feature these 
processes. In the two Higgs doublet version of quark-lepton symmetric models, tree-level 
flavour changing processes are, by contrast, unavoidable. As the lagrangian is symmetric 
between quarks and leptons, the flavour- changing leptonic terms are dependent on quark 
masses, and vice-versa, so leptonic processes are of more interest as the heavy quark masses 
give rise to proportionally large branching ratios. 

There are in fact a number of alternative ways of incorporating two Higgs doublets into 
quark-lepton symmetric models, but here we have concentrated on the model which gives 
the larger contribution to flavour- nonconserving decays, the Model 1 of Ref. 0, the main 
results concerning which will be briefly stated in Sec. II of this paper. Other models, such 
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as the Model 2 of Ref. ||, lead to leptonic Yukawa couplings with a dependence on lepton 
masses as well as quark masses. As stated above, it is the couplings proportional to quark 
masses that are of interest, as these lead to flavour- changing reactions. 

Thus the main focus of this paper will be on flavour- changing leptonic decays, in par- 
ticular those of muons since these give the most opportunity for experimental investigation. 
First in Sec. Ill the neutral Higgs mediated decay \T — > e~e + e~ will be investigated at tree- 
level, as will the charged-Higgs contribution to the decay pT — > u^e~i7 e . In Sec. IV a one-loop 
calculation of the process /i — > will be performed. This will lead on to an examination of 
the effect on the muon anomalous magnetic moment. Finally, the effect of the quark- lepton 
symmetric 2HDM contribution to b — > s'j will briefly be considered. 

II. YUKAWA COUPLINGS IN A QUARK-LEPTON SYMMETRIC MODEL 

The symmetry-breaking sector in the two Higgs doublet version of the quark-lepton 
symmetric model we will consider here consists of Higgs fields with the following quantum 
numbers under the symmetry group G q g @: 

Xi ~ (1, 3, l)(-2/3), X 2 ~ (3, 1, l)(2/3), 0! ~ (1, 1, 2)(1), 2 ~ (1, 1, 2)(-l), (1) 

where Xi X2 and 0i <-> 02 under the discrete quark-lepton symmetry. 

Symmetry breaking occurs in two stages. In the first stage, X2 acquires a non-zero 
vacuum expectation value (VEV), which breaks the leptonic colour group SU(3)£ down to 
SU(2)', and thus also breaks the discrete symmetry. The many other interesting features of 
quark-lepton symmetric models, such as the exotic leptons of the unbroken SU(2)' group, 
and the new gauge bosons which appear, are dealt with elsewhere 0. To achieve electroweak 
symmetry breaking, both 0i and 02 acquire VEVs. The corresponding Yukawa Lagrangian 
is 

£ Yuk = H(Ql) c QlXi + {FlYFlXi] + A 2 [M%?Xi + M^rX 2 ] 
+A 1 [Qld R <j) 1 + FEv R <j) 2 ] + A'JQlc^ + FEv R <Pl] 
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+HQlu r <P c 1 + F L e R <p c 2 ] + A' 2 [QlUr4> 2 + F L e R (f)i) + H.c, (2) 

where Ql and F L represent the left-handed quark and lepton weak doublets respectively, 
and cf) c is the charge-conjugate of the scalar field, <fr c = ir 2 <j)*. 

Electroweak symmetry breaking gives rise to physical fields from the two electroweak 
doublets as in normal 2HDMs ||. Thus there is a charged field H ± , a CP-odd neutral field 
7], and two CP-even neutral fields, h\ and h 2 . There is also one neutral physical Higgs field: 
that part of \ 2 which survives after quark- lepton symmetry breaking. Following Ref. [0, 
we will neglect mixing between this field and 77, h\ and h 2 , since \2 is expected to be much 
more massive than any of these electroweak fields. 

The resulting Yukawa interactions for these physical fields are then || 

£y u k = -d L m e u R H- + -u L m® irac d R H + + -u L m u e R H + + -e L m d u R H' + H.c, (3) 
u u u u 

% % % — % 

^Yuk = —7=-u L m e u R r] + —7=-e L m u e R r] + —7=-d L m^ 1Tac d R r] + —7=-v L m d v R r} + H.c, (4) 
y2u y2u y2u v2u 

^Yuk = —T^T^l m u {cos(uj - (p)fn + sin(cu - tp)h 2 } 



V2u 



+m e { — sin(u — ip)hi + cos(u — ip)h 2 } 

1 



u R 



m u {sm(uj — tp)hi — cos{uo — <p)h 2 } 



+m e {cos(cG> — (p)hi + sin(u; — (f)h 2 } 

1= d~h md{cos(u — <p)h\ + sin(u; — ^>)h 2 } 
\/2u L 

+m° irac {sin(u; - ip)h\ - cos(u; - tp)h 2 } d R (5) 
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rrid{— sm(u — <p)h\ + cos(a; — p)h 2 } 



y/2u 

+m° irac {cos(co> - y)h\ + sin(u; - (f)h 2 } v R + H.c. 



where u = \Ju\ + u 2 , with U\ and u 2 being the VEV's of <pi and <p 2 , respectively. Also in 
the above, tana; = W2/W1, and (p is a mixing angle relating the mass-eigenstate fields hi and 
h 2 to the CP-even parts of 0i and <p 2 . 

These Yukawa lagrangians are written in terms of the gauge eigenstate fermion fields, 
eg. u for the charge 2/3 quarks. To obtain the corresponding mass-eigenstate fields, U, the 
unitary diagonalisation matrices V£ R are introduced: 



U LjR = V^ R u LjR (6) 

and the diagonal mass matrix M u = V^m u V R . Analogous relations of course apply to the 
other kinds of fermions, d, e and v. 

Writing the Yukawa lagrangians in terms of the mass-eigenstate fields [and using the 
notation V[ 9 R = Vl R V£ R , so that for instance the Cabibbo-Kobayashi-Maskawa (CKM) 
quark mixing matrix is V£ ], we obtain for the flavour-changing interactions 

-U(yZ u M„V R iut 1+ + V R ie M e Vth-)DH + 

+^N(yr ] M u V£ e 1+ + V^M d Vtl-)EH+ + H.c, (7) 
-^-U(vrM e V^ 1+ )Ur, + -^-E(vr'M u V^ + )E v 
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Yuk 



''-Yuk 



rh 

^Yuk 



-^-D^M^^Dr] + -±-N(v^M d V^ + )N V + H.c. 

^u(vr e M e v^ 1+ )u 



(8) 



y/2- 



-L-E(vr j M u V^ 1+ )E 



sm(u — (p)hi + cos(u; — Lp)h 2 
sin(u — <f)hi — cos(ciJ — ip)h 2 



V2 



u 



D{yt v M v V R ^+)D sin(cu - if)ht - cos(w - tp)h 2 



(9) 



cos(a; - (p)h 2 

\/2u v / l 



+ H.c. 



It is important to observe that down-quark sector neutral flavour-changing vertices are 
proportional to neutrino Dirac masses. Because of the severe experimental upper bounds 
on neutrino masses, and because they are Dirac particles in our model, down-quark sector 
effects are negligible. This is pertinent because the constraint from K° — K° mixing would 
otherwise be very stringent (see Ref. for further discussion). We will thus concentrate on 
processes in the charged-lepton and up-quark sectors. 

In the rest of this paper, the following notation will be used: 

Oi,m 2 ,u 3 ) = (u,c,t); (d 1 ,d 2 ,d 3 ) = {d,s,b); 
(ei, e 2 , e 3 ) = (e, r); (u u u 2 , u 3 ) = (z/ e , z/ r ); 



m 



Vf 



2 -; z f 

H 



m 



f (i>2) 

2' W f 



m. 



2 ' 
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1 _ sin 2 (uo — if) cos 2 (cj — </?) 



m± 



+ 



m 



hi 



m 



hi 



mi 
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wf = 



mj_ 
m± ' 
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III. LEPTONIC DECAYS AT TREE LEVEL 



A. \T 



e e 1 e 



At tree level, the flavour- changing decay fi — > e e + e can be mediated by the neutral 
Higgs particles, rj, hi and /i 2 . The branching ratio for this decay is (with the approximation 



m 



<m 2 ), 



4 L 4m_ 
1 



fc,Z 



Em m V ue *V ue V ue *V ue 
m Uk m Ul v Rkl v Lk2 v Rll v Lll 

k,l 



(\S^m m V ue *V ue V ue *V ue 2 -\- \ m m V ue *V ue V ue *V ue 2 \ 
\ \ m u k ' n u l VRkl v Lk2 v Lll v Rll + | m u h m u l VLkl v Rk2 v Rll v Lll ) 

k ,/ k . / 

"I I 2 2~ m iifc m !ii m !im m «n I ^Lfcl ^Rfc2 *\L«1 ^ffll ".Rml ^Lm2 ^Lnl "iinl 



t 1 t k,Lm,n 



ii/ue*j/ue i/ue*i/ue i/ue* j/ue i/ue*i/ue , \\ 
"+" ^iifel K Lfc2 ^ffll ^Lil "Lml "flm2 ^ftil V Lnl "+" CX J J 

1 ?7le / 1 



+C.C 



))] 
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The above branching ratio involves many unknown parameters in the Higgs boson masses 
and the mixing matrices. In order to get an indication of how large B(/j, — > eee) might be, 
we will suppose that all of the mixing matrices display a hierarchical structure similar to 
that of the CKM matrix, which can be written in the qualitative form 



lee 6 



V 



e 1 e 2 
e 3 e 2 1 

where for the CKM matrix for instance, e ~ 0.22. The branching ratio is then approximately 



V 



(14) 
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eee 



m yue*yue t yue*yue , yue*yue \ 



m: 



+ 6 



m yue*yue I yue*yue , yue*yue * 

„ rn c V R21 V L22\" l uV Lll v Rll + '' t c» / L21 ^21, 



+ 



+ 



R 



R 



m M m_m 
(L i?) 



i fro 2 T/ ue *l/ ue T/ ue *T/ ue fro V ue *V ue -4-ro V ue *V 



R2l) 2 + CX 



(15) 



This branching ratio has the experimental upper bound B(fi — ► eee) < 1.0 x 10~ 12 . 

At this point, in order to get some idea of the relation of the branching ratio to the Higgs 
boson masses and mixing matrix elements, it is useful to assume that the matrices are of the 
form Eq. ([14]), and, for convenience, that they share the same e parameter. Of course, this 
loses some of the features of the unitarity of the mixing matrices, as well as the possibility 
of accidental cancellation, but nevertheless at least a qualitative understanding should be 
attainable. Eq. (|1|) can then be written: 



2\2 



B(fM — > eee) ~ ^ m c e ( m u + m c£ ) 



1 



2m 4 _ 



+ 



3mi 



+ 



m M m?.m^ 



(16) 



There are three extreme possibilities for the relations between the masses of the various 
neutral scalars, each of which leads to a slightly different relation between the mass and the 
e parameter. 
If m v ps m hl ps m h2 , 



If m 2 < m 2 hl ,m 2 h2 , 



If m 2 > m 2 hi ,m 2 h2 , 



— < — 

mi mi 



mj 



(17) 



m: 



m: 



m; 



(18) 



m: 



m5 



m 



2 ' 



(19) 



The dependence of the branching ratio on e for various choices of the Higgs boson masses is 
shown for these three cases in Fig. ([!]). The minimum value for the mass used, 48 GeV, is 



that obtained by the ALEPH collaboration ||. Only two lines for each e value are visible, 
because the plots for the conditions of Eq. ( |T8| ) and Eq. (|1^) are inseparable on the scale 
used. It can be seen from this that an increase in the precision of the measurement of the 
branching ratio by only a couple of orders of magnitude opens up for investigation vast new 
regions of the parameter space. Indeed, a considerable range is already excluded. It can be 
seen that for e = 0.22, corresponding to the CKM matrix, the Higgs boson mass has to be 
greater than about 330 GeV. 

Similar calculations can also be performed for the tauon as well, in such processes as 
t~ — > yU^/i+yU - and t~ — > e~e + e~. Writing their branching ratios in the same form as 
Eq. (|16D, we obtain where B T = B(t — > v T ev e ) = 0.1793 

B(r m) ~ Br ^K + ^f{^_ + ^4 + ^4 ( ^ + 3^ ' m 
B(r ee-e) ~ B r \ m ^(m c + m t ^_ + ^ + (21, 

Using Higgs boson masses of the order ~ 100 GeV, and for the e parameter e ~ 0.1, these 
branching ratios are of the respective orders 10~ 8 and 10~ 14 (here, as elsewhere in this 
paper, the value of m t = 150 GeV is used for the top quark mass, as this is in agreement 
with current experimental data [[]]). Experimentally, these branching ratios are less than 
1.7 x 10 -5 and 2.7 x 10~ 5 respectively, therefore these processes do not constrain our model 
to nearly the same extent as does \i~ — > e~e + e~ . 



B. n -¥ v^e v e 

At tree level in the Standard Model, the decay \i~ z/ A1 e~z7 e is mediated by the W boson. 
In our model it can also be mediated by a charged scalar, H ± . This can give a constraint on 
the various parameters (ie. mass and mixing matrices), since the experimentally measured 
muon decay agrees so well with Standard Model predictions. That is, the charged Higgs 
contribution must be of the order of the uncertainty in the decay rate, T(fi —>■ z/ M ez7 e ) = 
(2.99592 ±0.00005) x 10- 16 MeV §. 
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Again with the approximation m 2 e <C m' \, and using hierarchical mixing matrices as in 
Eq. (|T5D, the total decay rate, including both W and Higgs contributions is 



r — rv + + ^hw 



1 + \{yc\V^\ 2 + y t \V^\ 2 ) (y c \V^ 2 \ 2 + y t \V^ 2 f) 

x (mcyssyssyzs, + mtV&{v&v%> + y£jyz&))] , (22) 



where 



w 



(23) 



1927T 2 ' 

In Eq. (p2|), is essentially the leptonic equivalent of the CKM matrix, and must be used 
in the W boson couplings. 

Putting typical values for all the parameters (to# = 100 GeV, e = 0.1), the total contri- 
bution from the terms Y H and T HW is (T H + Y HW )/Y W = 1.3 x 10 -8 , which is over three 
orders of magnitude less than the experimental uncertainty, AT/T\v ~ 1-8 x 10~ 5 . 

At this point it should be noted that the minimum mass for the charged Higgs boson used 
in this paper is 45.3 GeV, as obtained by the ALEPH collaboration |J with the assumption 
that the branching ratio B(H + — ► t + u) = 1. In this particular 2HDM, with its quark-lepton 
symmetry, leptonic decays are proportional to the square of quark masses (in particular the 
top quark mass for t + u), whereas hadronic decays are proportional to the squares of lepton 
masses. Thus here we can take B(H + — > t + u) = 1, and so use the limit mjj > 45.3 GeV. 



IV. ONE-LOOP PROCESSES 

A. fi — ► e7 

One-loop processes such as fi — > and b — > have been extensively studied, both in 
terms of the Standard Model @ and 2HDMs ||. In the two Higgs doublet quark-lepton 
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symmetric model, the branching ratio for the decay /i — > is approximately [from Eq. (|A1|) 
of the Appendix, using the fact that neutrino masses are very small, and the approximation 
m\ < ml}: 



m. 



7/)- ?/i-T/ 1ie *T/ Me v ue *V ue 
w u k w u l v Lkl v Rki v Lli V R12 



11 v ue *V ue A -?/_, v de *V de — ?/)+ v ue *v 

III,, 



ue 
Rk2 



m a 1 



(24) 



Once again, the hierarchical structure Eq. ([TJ|) can be used. If the neutral scalars are all 
much heavier than the charged scalars, then the branching ratio Eq. (|24]) can be written 



£(//-> ery) 



a 



967T 



" l c v R21 v R22 + '' t * V R31 v i 



t K ii31 v R32 



(25) 



Alternatively, the neutral bosons could be much lighter — then the same possibilities as 
used for B([i — > eee), Eqs. ( ff7L |TS| , [T9|), can be used again: 
a 1 



#0 -> ery) 
B(j*-> ery) 

B(j*-> ery) 



24vr ml 
a 1 r 



m 2 V ue *V ue 4- m 2 l/ ue *l/ Me 
m c K L21 V L22 + m i ^LSl |/ L32 



+ 



i? 



967T 



2 (] rue*] rue q] rue*] rue t rue*] rue \ 
rn c\ V L21 v L22 vv L21 v R22 v L22 v R22J 



m 



j_ rrn 2(]rue*]rue q t~] rue*] rue ] rue*] rue s 
' t I ^LSl V L32 y IT~ V ^31 v fl33 V L33 v i?32, 



+ 



R 



a 1 



967T 



2 /True*] rue i Qi/ae*i/ue ]rue*]rue \ 
" l c\ V L2l v L22 yv L2l v R22 v L22 v R22/ 



m 



_\ ri - Y ~ i 2 / t rue* t rue . q T ] rue*] rue t rue*-\ rue N 

"" t ~''M v L31 ,/ L32 "r y "~ v L31 v fl33 V L33 V ,R32, 



+ 



(26) 



(27) 



(2f 



Experimentally, _B(/i — > e7) < 4.9 x 10~ n ||. Parameterising as per Eq. (0), the 
branching ratio has been plotted in Fig. (|2]) as a function of e for various values of the Higgs 
boson mass for the case in which the charged Higgs boson contribution is dominant. The 
corresponding plots for the other cases are similar. In all these plots the experimental limit 
is only a couple of orders of magnitude greater than the branching ratios corresponding to 
that region of parameter space of greatest interest. As for current limits, an assumed value 
of e = 0.22 gives the bound m# > 92 GeV on the mass of the charged Higgs boson. 
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Again, these calculations can be performed for the tauon as well as the muon. Once 
more using the e parametrisation, we obtain the branching ratios 



a 
96tt' 



,4,4 



^m r 



f 1 1 



4 



1 m,, 1 1 \ 2 
- h 9 

2 2 1 2 

mi m T mi mi 



r^J_ + 9 J_ xa - 



mi 



a 



cy) ~ -B T mfe e 

96vr * 



mi 



m 



m T m 

9 x2 
~~2 



m_ 



in 



VTLa 



m 



9 n2 
~2 



(29) 
(30) 



These branching ratios are ~ 10 8 and ~ 10 10 respectively (with masses ~ 100 GeV and 
e ~ 0.1), compared to the experimental upper bounds 5.5 x 10 -4 and 2.0 x 10 -4 . 



B. Muon Anomalous Magnetic Moment 

The calculation of the anomalous magnetic moment of the muon to one-loop is very 
similar to the above calculation of \i — > ej. The total contribution of the various Higgs 
particles is 

-9 £ — v^« Re (^^^r^f 2 )] (31) 

If typical values of the masses and e are substituted into Eq. fl3"ID , the anomalous magnetic 
moment is found to be less than ~ 10~ n . Experimentally it is known || that —13 x 10~ 9 < 
Aa„ < 21 x 10~ 9 , so this constraint is easily satisfied. 



C. s 7 

The nonleptonic decay b —>■ sj is currently of interest due to the recent detection of this 
process by the CLEO collaboration |9[], and is as yet not in disagreement with Standard 
Model predictions. The 2HDM contribution has previously been the focus of much investi- 
gation !|, and so it is worthwhile considering how the quark-lepton symmetric two Higgs 
doublet model differs in its effect on this process. 

11 



Using the notation of Grinstein, Springer and Wise |10| , the form factors C 7 (mw) and 
Cg(mw), referring to the operators for the photonic and gluonic vertices respectively, at the 
mass-scale Mw, need to be calculated. 

In the Standard Model, these form factors are given by 



C 7 (m w ) 
C 8 (m w ) 



~D(xt) 



where Xf = rrit/m^r, and 



A(x) = x 



D(x) = \ 



2l -v. 2 I 5 y, _7_ 

3 _ 12 12 

(x - l) 3 



|x 2 — x)lnx 



}_ rp 2 5 

2 2 



1 



+ 



(x - l) 4 
3x Inx 



(32) 
(33) 

(34) 
(35) 



(x-1) 3 (x-1) 4 - 
In this quark-lepton symmetric two-Higgs doublet model, the dominant contribution 
comes from the heavier third generation, in both the quark and the lepton sectors. Perform- 
ing the necessary one-loop calculation, again similar to that for \x — > ej, the result is (for 
simplicity only considering charged Higgs) 



1 1 7/ V de V de * 

C,( mw ) = -\a( X <) - i^^M« 

A Vt VLZ2 

1 1 ii v de v de * 

C s (m w ) = - l -D(x t ) - \ V -^^D{y t ) 



(36) 
(37) 



2 — 6y t V£g 

Because my/ = 83 GeV, mu > 45.3 GeV, and A(x) increases with x (ie. decreases with 
m\y or ma), A(y t ) < 3.36A(x t ). If the mixing matrices are similar to the CKM matrix then 
^L23^L33/^L32* ~ 1 ( & t an y ra te it cannot be too much greater than 1). The deciding factor 
is then y T /yt = m^/mf < 1, so the charged Higgs contribution should be negligible to that 
from the Standard Model. As stated above, this is in agreement with current experimental 
data. 



V. CONCLUSION 



If the two Higgs doublet version of the quark-lepton symmetric model is correct, then as 
we have seen, its contribution to various leptonic decays should be measurable for a wide 
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range of parameters, with only a small improvement in detector precision. As however the 
Higgs boson masses and the mixing between lepton families are unknown, it is possible that 
extreme values (i.e. very diagonal mixing matrices, or very heavy scalars), could lead to 
extremely small branching ratios. Measurement of the above decays is at least one way to 
obtain information about these parameters. 
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APPENDIX: 

The total decay rate for the process fi — > to one loop, including both Standard Model 
and charged and neutral Higgs contributions, is (assuming small lepton masses) 



-L /I lib,, ■ l j 



de 
L12 



-\rn^y/w- k w-V^V^V^V^ 2 ) 



i,k,l 

+rnJ B ^y^yt k {V^V^V^ + m^^yiy^V^V^V^ 



-m^wiw+y^y^rv^ - -m e r B ^wi k wiyr k lv^y^:v^ 2 
i 



--m e j^w- k w-y^iv^y^:v^ 



where Fo is the Standard Model contribution 



(Al) 



F 2 « - X -f (A2) 
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with Xf = mj/m^,. The other terms in Eq. ( |A1|) are 

rl rl—u ] 

I a = / rfw / cfo (1 — u — v)v/D w — , for small lepton masses, 
Jo Jo 12 

/•l rl—U | 

Ib = du dv (1 — u — v)u/D w — , 
Jo Jo 12 

i-l rl—u | 

h= l du I dv (\ — u — v) I D & —, 
Jo Jo 2 

1 /-l-u 



/' = / du / U du(i - u -v)v/D' w -, 
jo jo 6 

/•l rl—U j 

I' B = du dv (1 — it — v)u/D' ps -, 
Jo Jo 6 

1'=/ / dviu + v)/^ 

Jo Jo 2 

D = (1 - u - f - 1 - - ww e ) + 1, 

7J' = (1 - M - f )(1 - U>e. - DW,, - MW e ) + W7 6i . 
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FIGURES 

FIG. 1. Plot of the branching ratio for the process fx — > eee against the mixing matrix parameter 
e, for a range of masses of the neutral Higgs particles. The upper line in each pair corresponds to 
the mass m^, and the lower line to m rj and m^. The experimentally obtained upper limit on the 
branching ratio is also indicated. 

FIG. 2. Plot of the branching ratio B(n — > ej) against e for various values of the mass of the 
charged Higgs boson H + , with the current experimental upper bound indicated. 
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